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ABSTRACT 

Spitzer images resolve the debris disk around 7 Ophiuchi at both 24 and 70 /im. The resolved images suggest a disk radius 
of ~520 AU at 70 /im and >260 AU at 24 /im. The images, along with a consistent fit to the spectral energy distribution of 
the disk from 20 to 350 /im, show that the primary disk structure is inclined by ^50° from the plane of the sky at a position 
angle of 55°±2°. Among a group of twelve debris disks that have similar host star spectral types, ages and infrared fractional 
luminosities, the observed sizes in the infrared and color temperatures indicate that evolution of the debris disks is influenced by 
multiple parameters in addition to the proto-planetary disk initial mass. 

Subject headings: circumstellar matter - infrared: stars - planetary systems - stars: individual (7 Oph) 



1. INTRODUCTION 

Planetary debris disks are one of the best means to ex- 
plore the evolution of planetary systems. The Spitzer mission 
has made fundamental contributions to our understanding of 
them, including documenting the wide variety in properties 
such as amounts of excess emission, disk size and grain prop- 
erties derived from mineralogical features (Rieke et al. 2005, 
Su et al. 2006, TrilHng et al. 2008, Chen et al. 2006). Two 
possibilities for this variety are: 1.) disks share a similar evo- 
lution, but with a wide range of initial masses (e.g., Wyatt et 
al. 2007b; Kenyon & Bromley 2008); and 2.) there are sub- 
stantial evolutionary differences such as large collisions that 
dominate the amount of debris for a period of time (Wyatt et 
al. 2007a). There may also be a mixture of these possibilities 
(Rieke et al. 2005; Wyatt et al. 2007b). In any case, under- 
standing debris disks requires determination of their range of 
behavior where critical boundary conditions that might influ- 
ence these systems are held fixed. 

Here we report a detailed study of the debris disk around 
7 Ophiuchi, a 184 Myr-old AO V star (Song et al. 2001, Rieke 
et al. 2005) at a distance of 29.1 ±0.8 pc with an infrared ex- 
cess detected in the IRAS point-source catalog (Sadakane & 
Nishida 1986). Fajardo-Acosta et al. (1997) report marginal 
evidence that 7 Oph is extended at 60 iJ,m. In this paper, 
we describe imaging with Spitzer that confirms this result by 
clearly resolving the disk at both 24 and 70 /im. Although the 
radii of other resolved disks (e.g., Fomalhaut) are typically 
only about 150 AU, we find that the 7 Oph disk is almost 
four times as large. We combine this size measurement with 
the detailed spectral energy distribution (SED) of the disk to 
derive its properties, and compare them to those of a small 
sample of debris disks with similar host stars. 

2. OBSERVATIONS AND DATA REDUCTION 

The observations presented here are from various Spitzer 
programs. The IRS spectrum of 7 Oph was previously pub- 
lished based on an older pipeline reduction (Sll) by Chen 
et al. (2006). Here we present the S15 reduction with up- 
dated flux non-linearity corrections. Observations at 24 /im 
were obtained at two epochs in standard small-field photome- 
try mode using 3 secx4 cycles at 5 sub-pixel-offset positions, 
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resulting in a total integration of 840 sec on source for each 
epoch. The observation at 70 /im was done in fine-scale mode 
with 3 sec x 1 cycle at 4 sub-pixel-offset positions, resulting 
in a total integration of ^100 sec on source. A MIPS SED- 
mode observation was obtained with 10 sec x 10 cycles, and a 
1' chop distance for background subtraction. 

All of the MIPS data were processed using the Data Analy- 
sis Tool (Gordon et al. 2005) for basic reduction (e.g., dark 
subtraction, flat fielding/illumination corrections), with ad- 
ditional processing to minimize artifacts (see Engelbracht et 
al. 2007). For the 70 /im fine-scale data, each of the on-source 
exposures was subtracted from an average background com- 
posed of before and after off-source exposures to remove the 
time-dependent column offsets in the background. The final 
mosaics were combined with pixels half the size of the phys- 
ical pixel scale (Figure [TJ. The MIPS SED-mode data were 
reduced and calibrated as described by Lu et al. (2008). 

At 24 /im, the non-color-corrected measurement gives 
434±3 mJy using an aperture radius of 15" with a background 
annulus of 30"^2", and an aperture correction of 1 . 143. The 
expected photospheric flux density in the 24 /im band is 236 
mJy, resulting in a 24 /im infrared excess of 202 mJy after a 
color correction of 1.021 (for a blackbody temperature of 80 
K, which is the characteristic color temperature of the disk 
(see §3.2)). The total color-corrected flux density in the 24 
/im band is 438±18 mJy, including the uncertainty in the cal- 
ibration factor of 0.0454 MJy sr"' (DN s"')"' (Engelbracht et 
al. 2007). For the 70 /im fine-mode photometry, we use an 
aperture with a radius of 16", sky annulus of 18"-39", and an 
aperture correction of 1.9. A total flux density of 1140 mJy 
was estimated using a flux calibration factor of 2894 MJy sr"' 
MIPS70F-' (Gordon et al. 2007). The expected photospheric 
flux density in the 70 /im band is 26 mJy; therefore, the disk 
flux density is 1 185 mJy after a color correction of 1.064 (for 
a blackbody temperature of 80K). The total color-corrected 
flux density in the 70 /im band is 1211±121 mJy. All avail- 
able photometry and spectra of 7 Oph are shown in the SED 
plot in Figure |2] 

3. ANALYSIS 

3.1. PSF Subtraction 

The excess emission at 24 /im is roughly equal in bright- 
ness to the expected photospheric value and yet the source 
appears point-like at 24 /im. To reveal the disk, the Point 
Spread Function (PSF) must be subtracted. We used deep ob- 
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Fig. 2. — Spectral Energy Distribution (SED) of 7 Oph. The photometry 
of 80±20 mJy at 350 /im is obtained with the Caltech Submillimeter Ob- 
servatory (Dowell et al. in prep.)- The optical and near-infrared photometry 
is from the Simbad database and 2MASS catalog, while the other symbols 
are defined on the plot. The dotted-dashed line is the output SED using the 
best-fit parameters. 
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Fig. 3. — The excess SED of the 7 Oph disk. Symbols and lines are the 
same as in Figure |2] except for the stars showing the ten photometric points 
used for fitting. Two single component models are shown (p=0 with grains 
of 5-63 fim and x^=0.8, and p=l with grains 1.6-100 /im and x^=10.4). As 
an illustration, the total SED of a three-component (not unique) disk model 
with x^=0.7 is shown along with the SED of each component. 



servations of the non-excess A4 V star T3 Eri to determine the 
PSF. At 24 /im, the Full- Width-Half-Maximum (FWHM24) of 
the T3 Eri image is 5. "52x5. "48 using a 2-D Gaussian fitting 
function on a field of 26"with an uncertainty of 0"02x0"02 
based on ^^100 calibration star observations. The T3 Eri image 
was scaled to the expected photospheric flux density of 7 Oph 
and to 1.6 times this value for PSF subtraction. The higher 
value (over-subtraction) allows for the possibility of a signif- 
icant portion of the excess emission lying close to the star as 
in Fomalhaut (Stapelfeldt et al. 2004). The results of the PSF 
subtractions are shown in Figure [T] An elliptical outer disk 
along the position angle of ~56° is seen in the over-subtracted 
image, while the photospheric-subtracted image confirms that 
much of the disk emission lies within the PSF. The disk is at 
least 9"(--260 AU) in radius at 24 ^m. 

To test for any color dependence in this process, we sub- 
tracted the image of T3 Eri from the image of Lep, a de- 
bris disk known to be largely confined to within 3 AU (0."14) 
radius at 20 fim (Moerchen et al. 2007), with various scal- 
ing factors. There are no disk-like residuals in the subtracted 
image. The FWHM24 of the photosphere-subtracted C, Lep 
image (5."61x5"55) gives a characteristic resolution for a 
red point source, slightly broader than given by a blue PSF 
(r3 Eri). The FWHM24 of 7 Oph before photospheric subtrac- 
tion is 5. "84x5. "61, which is only 1.06 x 1.03 broader than a 
blue PSF. However, this difference is significant at 16-<t com- 
pared to a blue PSF and ll-cr to a red PSF The FWHM24 of 
the 7 Oph disk after photospheric subtraction is 6"41 x5."75, 
which is 1 . 14 X 1 .03 broader than a red PSF represented by the 
unresolved core of the C, Lep disk. 

The source appears slightly elongated at 70 fxm (Figure Ic). 
No photospheric subtraction is needed because the disk is 
~50 times brighter than the stellar photosphere. However, the 
faint extended disk emission is more evident after the subtrac- 
tion of a PSF representing the bright inner disk (Figure Id). 
The faint extended emission can be traced up to ~18"('^520 
AU, at 3-(T level) in radius along the position angle of ~55°. 
The disk has a FWHM70F of 18"5x 16."1 determined using a 
Gaussian fitting on a field of 55". Examining two blue cali- 
bration stars (HD 48348 and HD 48915), a typical FWHMvof 

MIPS 24 fim measurement (247±3 mJy) and an updated IRS spectrum 
show no excess in this source. 



of 15."3xl4"4 with an uncertainty of 0."2x0."06 was deter- 
mined, suggesting that the 7 Oph image core at 70 fxm is 
1.2x 1.1 times wider than a blue PSF, and resolved at 16- and 
28-0- levels for the major and minor axises. The position an- 
gle of the disk determined from both 24 and 70 /im images is 
55°±2°. The MIPS sHt (width of 19"8) was at the position 
angle of ~95° when the SED-mode observation was taken. A 
small fraction of the faint extended disk emission was outside 
the slit, which may explain the slightly lower (~5%) SED- 
mode flux compared to the fine-mode photometry. 

3.2. Spectral Energy Distribution Fitting 

The broadband excess photometry at 24 and 70 /im indi- 
cates a color temperature (T^.) of 8 1 K, which over predicts the 
observed flux at 28-35 /im and 55-65 /im (see Figure|2]). This 
suggests that the dust grains in the system have a wider range 
of temperatures. To model the disk SED (Figure [3]), we use 
a simple geometrically-thin debris disk model where the cen- 
tral star is the only heating source and the dust is distributed 
radially between inner (R,„) and outer (Rom) radii according 
to a r~P power law for the surface number density. The grains 
in the disk follow a n{a) ^ a"* size distribution with a mini- 
mum radius of «„„>, and a maximum radius of a,„a.K, assumed 
to be the same throughout the disk. To minimize the free pa- 
rameters, we assume that the grains are astronomical silicates 
(Laor & Draine 1993) with a density of 2.5 g cm"^ and are in 
theoretical collisional equilibrium; i.e. q = -3.5. 

We determined the best-fit Kurucz model for the stellar pho- 
tosphere by fitting the optical to near-IR photometry (Su et 
al. 2006), and estimated the stellar temperature and luminos- 
ity to be 9750 K and 28.6 Lq. With these stellar parame- 
ters, grains smaller than ~5 /im are subject to radiation pres- 
sure blowout since (3, the ratio between radiation and gravi- 
tational forces (Burns et al. 1979), is > 0.5. To capture the 
slope of the IRS spectrum and MIPS SED-mode data, we se- 
lected two points (26.5 and 30 /im) in the IRS spectrum and 
four wavelengths (55.4, 65.6, 75.8, and 86.0 /im) in the MIPS 
SED-mode data for the global SED fitting. Combined with 
the broadband photometry, a total of ten data points shown 
as stars in Figure [3] are used to compute to determine the 
best-fit debris disk model. 

To fit the fairly flat slope between 55 and 95 /im and the 
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Fig. 4. — Comparison cuts along the disk major axis at 24 /im (upper panel) 
and 70 /xm (lower panel). The cuts at 24 [im are averaged over a width of 
5", and a width of lOf'5 at 70 iim. A dotted line in each panel indicates a 
nominal l-cr detection on the image. 



observed flux at 350 /im simultaneously, the maximum grain 
size has to be <100 /im. This might be an artifact of the 
assumption that the grains have the same size distribution 
throughout the disk (see ® . In addition, since there are no 
obvious silicate features, the majority of the grains in the disk 
are larger than a few microns. Therefore, we assume «„„„ to be 
the blowout size, ^5 [im, and fl,„av to be ^^63 /im, and that the 
disk has a constant surface density consistent with a Poynting- 
Robertson (P-R) drag dominated disk. With these grain prop- 
erties and density distribution, the SED is consistent with a 
1 -component disk with an inner radius of 13;I;3 ^ '^'^'^ outer 
radius of 430;!;3o"AU. However, a P-R drag dominated disk 
is not expected physically given the high disk density (frac- 
tional luminosity fd ^9 x 10"^) since collisions should destroy 
grains before they can spiral inward (Wyatt 2005). 

We also explored the possibility of fitting the SED with an 
outflow disk, in which the majority of the grains are smaller 
than the blowout size and the surface density follows a r"' 
power law (e.g., the Vega debris disk; Su et al. 2005). Various 
combinations of a„,ar, Rm and Ro,„ were tried, and none 
yield a correct spectral slope between 20-30 /im (one example 
is shown in Figure |3]l. In addition, because smaller grains 
have higher temperatures than larger ones at the same radius, 
models including grains smaller than the blowout size force 
R,„ to be a large distance (~100 AU=3."4), contradicting the 
observed compact core of the disk. This suggests that the 
7 Oph system cannot be an outflow disk entirely. 

3.3. Debris Disk Imaging Fitting 

Limited by the spatial resolutions, no multiple disk com- 
ponents can be distinguished in the observed images except 
for the extent of the disk. To further constrain the model, 
we constructed model images with the best-fit SED param- 



eters {p=0, a„i„ r^S /im, a,„ax ^^63 /xm, R,„=13;I;3 and 
R„„,=430!^^° AU) from the SED fitting. These model images 
were then convolved with appropriate PSFs to compare with 
the observed images based on the FWHMs and radial surface 
brightness distributions. 

For a geometrically thin disk model, the observed disk mor- 
phology only depends on how the disk is viewed, i.e., the in- 
clination angle of the disk with respect to the plane of the sky, 
/. We first explored the influence of the outer radius by fixing 
R,„=13 AU, the grain pai-ameters and /=90°. The FWHM24 
of the convolved model images with outer radii ranging from 
350 to 1000 AU shows no change, but a major-axis FWHM70F 
of the convolved model images varies from 17."0 to 20"4 with 
the minor-axis FWHMvof of 14."2. This again confirms that 
the 24 /im emission arises from warm grains that are not sen- 
sitive to the outer extent of the disk. The observed major- 
axis FWHM70F of the disk (18"5) impHes a R„,„ of -520 
AU. With the R,„, Ro,„ and grain parameters fixed, we deter- 
mined the inclination angle by fitting both major- and minor- 
axis FWHMs. Convolved model images were constructed for 
face-on (/=0°) to edge-on (/=90°) disks in increments of 5°. 
The major-axis FWHMs of the model images remain the same 
within measurement errors at both 24 and 70 /im. From 0° to 
90°, the minor-axis FWHM24 and FWHM70F changed from 
6"! to 5"6 and 17."6 to 14."2, respectively. We conclude that 
the observed FWHMs at 24 and 70 /im are consistent with a 
disk viewedat/=50°±5°. 

For a quantitative comparison, we show cuts along the disk 
major axis for both the observed and best-fit model images at 
24 and 70 /im in Figure HI Similar cuts made on point-source 
PSFs, scaled to match the peak values of the observations, 
are also shown. Our simple 1 -component, bound-grain debris 
disk model is in very good agreement with the observed disk 
surface brightness at 24 /im over a dynamic range of two or- 
ders of magnitude. At 70 /im, the observed signal-to-noise 
ratios are not as good as at 24 /im, but the model agrees to 
within the estimated uncertainties and approximates the data 
much better than a point source. 

4. DISCUSSION AND SUMMARY 

Based on the Spitzer images and SED modeling, the disk 
around 7 Oph has inner and outer boundaries of —13 and 
—520 AU, respectively, viewed at an inclination angle of 
—50°, assuming a 1 -component constant-surface (p=0) den- 
sity disk consisting of bound (5-63 /tm) astronomical silicate 
grains. The dust mass is — l.Ox 10"^ M^. The star's projected 
rotation velocity is 210 km s"' (Royer et al. 2007), suggest- 
ing an equatorial velocity of 256-300 km s"' with / ranging 
from 55°to 45°. Although these values are very close to the 
maximum permissible velocity (Gulliver et al. 1994), the es- 
timated disk inclination angle is within the permitted range of 
plausible stellar inclination angles. 

The broadband SED alone provides too little information to 
constrain properly the variables needed for even the simplest 
disk models, but a better understanding can be achieved by 
considering other constraints. The exact R,„ of the disk is 
sensitive to a,,,,,, used in the model. The fact that the excess 
emission starts at <15 /im suggests that dust exists as close 
as —10 AU, assuming blackbody radiators. Since the grains 
in a debris disk are generated by collisional cascades from 
large (— km-size) parent bodies, a,„a_t is usually assumed to be 
— 1000 /im. To make the single component model consistent 
with the faint extended disk detected at 70 /im and with the 
assumed grain size distribution, a,„a_t has to be reduced to —63 
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/xm. A more physically realistic model will involve multiple 
components while requiring the presence of very large grains 
in the disk. 

We used a three-component disk to derive a model con- 
sistent with the expected grain physics and the observed 
spectrum/SED. The three components are: (1) a narrow 
((SR'^-^O.IR) ring of parent bodies at ^80 AU where large 
grains are located; (2) an outflow disk (80-520 AU) with 
density cx r"' outside the narrow ring consisting of grains be- 
ing ejected by radiation pressure with sizes smaller than the 
blowout size; (3) a component with density oc r° interior to 
the narrow ring and extending to ~10 AU (e.g., as of grains 
slowly drifting inward due to P-R drag). Many combinations 
of parameters in this three-component disk model would yield 
acceptable fits to the observed data. One example is shown in 
Figure |3] but it is not unique without further constraints. 

Given the exceptionally large disk seen in the 7 Oph im- 
ages, one may wonder whether it is another Vega-like disk 
(an outflow disk observed at 24 and 70 fim extending from a 
ring of large parent bodies that dominate the emission at 850 
^m, Su et al. 2005). The failure to fit all of the available data 
with our simple 1 -component outflow model suggests that the 
contribution by grains escaping from the 7 Oph disk is weak, 
unlike for the Vega disk where such grains dominate the disk 
emission (^100% at 24 /im and r^75% at 70 pim). Persua- 
sive evidence of 7 Oph being another Vega-like disk would 
require the disk size at 350 fim to be substantially smaller 
than the size inferred in the infrared. However, the signal-to- 
noise of the 350 /im data is too low to show the extent of the 
emission. Observations at longer submillimeter or millimeter 
wavelengths and higher spatial resolution are needed to deter- 
mine the exact location and amount of the cold component. 

An attractive comparison sample of debris disks for 7 Oph 
is defined by systems with stellar ages between 150 and 400 
Myr, spectral types between AOV and A3V, and between 
10-5 and lO"'*, for which Su et al. (2006) find a total of 12 
stars^ including the resolved disks around Vega and Fomal- 
haut. The disk around Fomalhaut defines the expected tem- 
plate structure given our knowledge of the solar system. It 
consists of two separate dust belts: a Kuiper-Belt-analog ring 

' Besides 7 Oph, Vega and Fomalhaut: l Cen: ~250 K, Lep: Tc ~2()6 
K, a CrB; ~125 K, /3 UMa: ~116 K, 30 Mon: ~105 K, HD 38056: 
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with a radius of ^150 AU prominent at A >70 /im and a 
fair amount of zodiacal-analog warm dust located <20 AU 
from the star. Little emission originates from outside of the 
Kuiper Belt ring. In the Fomalhaut disk, the ratio between 
the central unresolved and the extended emission is 3-6 at 
24 fj,m (Stapelfeldt et al. 2004), while the ratio is --2-3 in 
the 7 Oph disk. The Kuiper-Belt-analog ring in the Fomal- 
haut system totally dominates the emission at 70 /im, but at 
the distance of 7 Oph ('-~^4x that of Fomalhaut), such a ring 
would be unresolved. 

Two of the twelve systems Lep and l Cen) have 24 to 
70 /im color temperatures > 200 K (Su et al. 2006). Ground- 
based 20 /im imaging of ( Lep shows it to be dominated by a 
very compact inner warm disk (Moerchen et al. 2007). Three 
more systems (30 Mon, /? UMa, a CrB) have color temper- 
atures of 100-150 K, indicating relatively large amounts of 
material in their inner zones. All of these systems are likely 
to have more debris-generating activity close to the star than 
the disks around Fomalhaut (color temperature of 70 K) or 
7 Oph. Of the remaining seven systems, at least two have 
anomalous activity in their outer zones: Vega with its large 
(-800 AU) halo, and 7 Oph with a disk --520 AU in radius. 

It appears that much of the overall debris disk behavior can 
be explained as the result of similar evolution from starting 
points that differ in the initial disk mass (Wyatt et al. 2007b). 
However, the large variety in behavior among the sample of 
twelve middle-aged disks demonstrates that one or more addi- 
tional parameters must also be at play. For instance, the pres- 
ence and number of massive planets, the recent collisional his- 
tory of the system (Kenyon & Bromley 2005, 2006), as well 
as the variety in the initial disk size may all play important 
roles in the wealth of disk behavior. 
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Fig . 1 . — MIPS 24 and 70 images of 7 Oph with colorbars sliowing the surface brightness levels between the 1 -u detection and the peak value of the image, 
(a) shows the deep 24 /im image of 7 Oph after photospheric PSF subtraction, (b) shows the over-subtracted image of 7 Oph at 24 /im, where an elliptical disk 
is evident, (c) shows the disk image at 70 /xm. (d) shows the over-subtracted image at 70 /xm after applying boxcar smoothing. The contours show the detection 
boundaries at 1-, 3- and 5-cr levels; \-a =0.126 mjy arosec"^. AH panels have the same orientation with N up and E to the left. A circle at the comer of each 
panel indicates the typical resolution. 



